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Abstract 

Discrete transformation for 3- waves problem is constructed in ex- 
plicit form. Generalization of this system on the matrix case in three 
dimensional space together with corresponding discrete transforma- 
tion is presented also. 







1 Introduction 



The problem of three waves in two dimensions arises in different form in 
many branches of the mathematical physics. Its application to problems of 
radiophysics and nonlinear optics reader can find in 0. In connection with 
the inverse scattering method it was investigated in and considered in 
details in numerous further papers. 

The goal of the present paper is generalize this system on the space of 
three dimension with simultaneous exchanging the unknown scalar functions 
on the operator valued ones. The last generalization allows to include into 
consideration the quantum region with the Heisenberg operators as unknown 
functions of the problem. 

It is necessary to mention that numerous different (by the form) discrete 
transformation were used up to now with respect to two components hier- 
archies of integrable systems, which are connected with so called Darboux- 
Toda, Lotke-Volterra and Heisenberg substitutions. In the present paper 
we come to the substitution which connect six independent functions, which 
corresponds to A 2 algebra but not to A\ as it was in the case of the two 
components systems. This substitution may be considered as the integrable 
mapping, connected six initial function with the six final ones. Substitutions 
of the present paper do not coincide with recently introduced in || so-called 
Ultra- Toda mappings. And the last comment. 

The method of this paper without any difficulties can be generalized on 
the case of n-th wave problem. In this case the number of independent 
variables of substitution will be (n x (n+ 1) which coincides with the number 
of positive and negative roots of A n algebra. 

The traditional way for obtaining the system of equations for 3-waves 
problem in (1 + 1) dimensions is the L — A pair formalism 
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where c, d four arbitrary numerical parameters (J2 c i — 0, XMi — 0) and 
(A, B, D, E, P, Q) are unknown functions of the problem. The system of 
equations for them has the form 

(d 2 - dx)P x - (02 - c^Pt + vEQ = 0, (d 3 - dx)Q x - (eg - ci)Qt - z/PA = 0, 

(da - d2)Ar - (c 3 - c 2 ) A + z/QP = 0, [d x - d 2 )B x - (d - c&)Bt - zMP = 0, 

(1) 

(di - d 3 )D x - (a - c 3 ) A + /yEP = 0, (da - d 3 )^ - (c 2 - c 3 )P t - ^DP = 0, 

where v = 3(cid 2 — c 2 di). Soliton like solution for the last system of equations 
may be easily found with the help of the technique of the paper [|[]. From 
these results it can be found the first steps of the discrete transformation 
with respect to which the last system is invariant. In present paper the 
reader have to consider the form of the discrete transformation as the lucky 
guess. 



2 Discrete transformation 

We can assume (but this is a direct corollary of the results of [Q) that the 
three "new" functions (Q, A, P), denoted by the bar symbols, connected with 
the old ones as following 

Q = -, A = --, P = - 
^ D' D D 

satisfy the system ([]]). Then from the first and second equations of the first 
column it is possible to determine E, B functions with the result 

E = —. 1 ~ {D X - (c 2 - Cl )BE) + 1 {E x + (c 2 - Cl )DP), 
(ci - c 3 ) D (c 2 - c 3 ) 

B = . 1 r -{D x - (c2 - a) BE) - . 1 A B X - (ca - c 2 )AD). 

And in a self-consistent way determine from the second and first equations 
of the second column D. We will not present here this not very simple 
expression, because in few lines below we will have observable expression 
for this value. Straightforward but tedious calculations show that the third 
equation of the first column is also satisfied (D 3 = —BE). 
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For further consideration it is more suitable to introduce three dependent 
variables ( £ + rj + a = 0) 

£ = (d 2 —di)t+(c2—c 1 )x, r] = (d 3 — d 2 )t+(c 3 — c 2 )x, a = (d 1 -d 3 )t+(c 1 -c 3 )x. 

In each pairs of variables (£, rf), (£, a), (77, a) the differentiation operators take 
the form 

Q 2 = fcfa)^ - ^ld t \ = \ ds-d a -d a 

d (d^g , Ss-ssLq^J \d v -d c d v 

Really the explicit form the generators of differentiation via (£, 77, a) variables 
will be not essential. Now the system (p]) looks much more attractive 

P x = -QE, A 2 = -BQ, Q 3 = —PA 

B 1 = -AD, E 2 = -DP, D 3 = -EB (2) 

In the last form the system is obviously invariant with respect to permu- 
tation of the indexes of differentiation with the simultaneous corresponding 
exchanging of unknown functions. The discrete transformation of the be- 
ginning of this section may be rewritten in more symmetrical form (we will 
denote it with the help of the symbol T 3 ) 

1 BE 
® = D' A= ~D' P= D> 

B = D(^) 2 , E = -D(^) 1 , ^ = DQ-{\nD\ 2 

By the permutation indexes (1,3) (together with corresponding exchanging 
of unknown functions) it is possible to obtain the T\ discrete transformation 
with respect to which the system (fj) is also invariant 

P = I, Q= ± E= _D 
B v B' £' 

D = B(^) 2 , A = -B(^) 3 , ^ = BP-(\nB) 2i3 
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And at last the discrete transformation T 2 has the form 

— 1 — D — P 
A =E< B= E< Q --E- 

D = -E(^) 1 , P = B(|)a, ^ = EA-(\nE) lt3 

In the form presented above substitutions Ti may be considered as a map- 
ping, connected six initial (unbar) functions with six final (bar) ones. From 
the other side each substitution may be considered as the infinite dimensional 
chain of equations. For instance the corresponding chain of equation in the 
case of T\ substitution has the form 

r>n+l r>n T~\n An 

F-r = -(lnB n ) 23 , D n+1 = B n (—) 2 , A n+1 = ~B n ( — ) 3 (3) 

pin An 
E n+l = Qn+1 = £L 

B n ' ^ B n 

In the first row we have the lattice like system connected 3 unknown functions 
(B, D, A) in each point of the lattice. The first chain for B functions is exactly 
well known two dimensional Toda lattice. 



3 Some properties of the discrete transfor- 
mations 

All constructed above discrete transformations are invertible. This means 
that unbar unknown function may presented in terms of the bar ones. For 
instance T 3 _1 looks as 

D = i. B=-i, E=Z, 
Q Q Q 

P=-Q(^h, A = Q(~)i, ^ = DQ-{\nQ\ 2 

It is not difficult to check by direct computation that discrete transfor- 
mations Tj are mutual commutative (T{Tj = TfTi) on the solutions of the 
system (0). 
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We present below corresponding calculations to prove that T{T 2 = T 2 T\ = 
T 3 . Indeed result of the action of Ti on some solution of the system @ is 
the following 

B v B B 
& = B(^) 2 , A' = -B(^) 3 , ^- = BP-(\nB) 2j3 
Action of the T 2 transformation on this solution leads to 
. <>i 1 B „r,-i T) 1 „ . -B . 



<5 21 = -^T = ^ ^ = -EH^h = -^(BilnDh-B^ = QD*-D(\nD) 12 

P 21 = E\^h = §, = (^W - ^(ln^) 13 = -2>(§)i 

The same calculation repeated in the back direction shows that PF 1 ' 2 = 
VK 2,1 = W 3 - the result of application of the T3 transformation to an initial 
solution W. 

Thus from each given initial solution Wq = (A, P, Q, E, B, D) of the sys- 
tem (0) it is possible to obtain the chain of solutions labeled by two natural 
numbers (Zi, l 2 , or (I3)) the number of application of the discrete transforma- 
tions (Ti, T 2 , T 3 ) to it (as it was shown above T{T 2 = T 2 Ti = T 3 ). 

The arising chain of equations with respect to (D, B, E) functions are 
exactly two dimensional Toda lattices. Their general solutions in the case 
of two fixed ends are well-known ||. As reader will be seen soon this fact 
allows to construct the many soliton solutions of the 3-wave problem in the 
most straightforward way. 



4 Resolving of discrete transformation chains 
4.1 Two identities of Jacobi 

We begin from the following obvious equalities for determinants of n-th order 
Det n (T n ) = D n F"- 1 a ) = D n ^(T n ^)(r - bT^a) = D n _ x {T n _ x )f 
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where T n _i is (n — 1) x (n — 1) matrix, a, b are (n — 1) dimensional column 
(row) vectors respectively and r scalar. 

By the same reason the following formula takes place 



D r , 
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7 m - ¥T n _ 2 a L Ti2 - b L T n _ 2 a z 
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&1 T H T 12 = D n - 2 {T n - 2 )D 2 I ,2 T -1 .1 _ i.2 T -l „2 

V 6 2 r 21 r 22 ' Vr 2 i - 6 l n _ 2 a r n - b l n _ 2 



where a 1 , b l are (n — 2) dimensional columns (rows) vectors, r^- components 
of 2-th dimensional matrix. It is obvious how relations of these types may 
be continued. 

Now using results above let us transform the following expression 

n ( T n -\ a 1 \ / T n _! a 2 \ n / T n _! a 2 \ / T n _x 
H 6 1 nJM 6 2 r 2 J M 6 1 rJM 6 2 

D n _ x {T n „ x )D 2 ^ tm _ h 2 T -^ a i Tu _ b 2 T -i ia 2 ) ~ Ah-i 

We will treated the last equality as the first Jacobi identity. By the same 
technique it is not difficult to show that the following equality takes place 
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This equality we will use many times in what follows and will call second 
Jacobi identity. These identities can be generalized on the case of arbitrary 
semi-simple group. Reader can find these results in ||. 

4.2 Concrete calculations 

Let us take initial solution in the form 

Q = A = P = 0, B = B(2), E = E(1). D 3 = -BE (4) 
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Application to this solution each of inverse transformations T~ is mean 
less via arising zeroes in denominators. The chain of equations under such 
boundary condition we will call as the chain with the fixed end from the left 
*from one side). 

The result of application to such initial solution l 3 times T 3 transformation 
looks as (for the checking of this fact only two Jacobi identities of the previous 
subsection are necessary) 
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where A n are minors of the n-th order of infinite dimensional matrix 
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and A^,A^ are the minors of I3 order in the matrices of which the last 
column (or row) is exchanged on the derivatives of the corresponding order 
on argument 1 of E function (on argument 2 of the B function in the second 
case). 

In what follows the following notations will be used. W l3,11 , (W l3,12 )- 
the result of application discrete transformation T la T l1 (T l3 T l2 ) to the cor- 
responding component of the 3- wave field. A Zs,Zl (A /s,/2 ) - determinant of 
^3 + h {h + h) orders, with the following structure of its determinant matrix. 
The first Z3 rows (columns) of it coincides with matrix of (§) and last l\, (I2) 
rows (columns) constructed from the derivatives of B, (E) functions with 
respect arguments 2, (1). 

The result of additional application of l\ times T\ transformation to the 
solution (|5|) looks as 



M 3 ,Zi-l 
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Q(h,h _ ( ]_y 3 +h-1 ^3-l,U £}(Wl _ ^_]_y 3 +«i^3+lA ^ 

A Wl A Wl 

We do not present the explicit form for components VKQs,^, which can be 
obtained without any difficulties from (0) by corresponding exchanging of 
the arguments and unknown functions. 



5 Many-soliton solution of the scalar 3- waves 
problem 

The system (0) allows the following reducing ( under additional assumption 
that all operators of differentiation are the real ones d a = 5*) 

P = B*, A = E*, Q = D* (8) 

In this case the system (g) is reduced to three equations 

B 1 = -DE*, E 2 = -DB*, D 3 = -BE (9) 

for three complex valued unknown functions (E, B, D). 

Now we would like to demonstrate how the multi-soliton solutions of 
the system (0) may be obtained with the help of the technique of discrete 
transformation in the most straightforward way 

With this aim let us consider the action of the direct and inverse T i; T~ x 
transformations on the reduced solution of the system (|9]). The trick consists 
in the fact that discrete transformation does not conserve the condition of 
the reality (§) and starting from the solution of the reduced system we come 
back to solution of irreductible one and in some cases vice versa. We will 
denote the three dimensional vector (Q, P, A) by the single symbol Q and the 
by symbol D three dimensional vector (D, B, E). Then the result of actions 
of direct and inverse transformations on solution satisfying the condition of 
reality Q = D* is the following 

TT&D*) = {ti) n {q,d), T~ n (D, D*) = (d*,q*) 
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where U are point like symmetries of the system (0) 

t 3 (Q, P, A, D, B, E) = (Q, -P, -A, D, -B, —E), 

t 2 (Q, P, A, D, B, E) = (-Q, -P, A, —D, -B, E), 

t 3 (Q, P, A, D, B, E) = (-Q, P, -A, —D, B, -E) 

It is obvious that tj = 1. Thus if we apply 2n times discrete transformation 
to initial bad (nonreduced) solution (0, D) and as a result obtain (t n 
vecD*, 0) then in the middle of the chain we will have solution satisfying the 
condition of reality, which coincides with n soliton solution of the reduced 
system @. 

— * 

The solution of the chain with the boundary conditions Q = on the left 
end of the chain and D = on the right side we will call as the chain with 
fixed ends. Really condition D = is the system of equation from which 
initial functions D, B, E (see (|D m ay be defined as the solutions of ordinary 
differential equations (see Appendix II). 

6 Matrix three waves problem in the space 
of three dimensions and its discrete trans- 
formation 

In all calculations above we have never used (except of concrete resolving of 
discrete transformation chains) the condition that operators of differentiation 
are connected by the condition 

di + d 2 + d 3 = 

as it follows from the definition of this operators. So we can consider the 
system (0) where all three operators are independent from each other and 
correspond to differentiation with respect to one of coordinates of three di- 
mensional space. The second generalization consists in possibility to consider 
the unknown function in (fj) as the operator valued ones. Of course in this 
case the order of the multiplications are essential and exactly coincides with 
fixed by the formula (0). 
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T3 discrete transformation in this case looks as 

Q = D~ 1 , A=-BD-\ P = D- 1 E, 

B = -D(BD- 1 ) 2 , E = —D(D~ 1 E) 1 , D~ l D = QD - (D~ 1 D 2 ) 1 
By the same technique for T\ we have 

P = B~ 1 , Q = B~ l A, E = -DB-\ 

D = (DB- 1 ) 2 B, A = -B(B- l A) 3 , BB^ 1 = BP — (B- i B~ 1 ) 2 
And at last the discrete transformation T 2 looks as 

A = E~\ B = E- 1 D, Q = -PE-\ 

D = -E(E- 1 D) 1 , P=(PE- 1 ) 3 E, E~ X E — AE — (E~ 1 E 3 ) 1 

As in the scalar case the discrete transformations in the case under con- 
sideration are mutually commutative. The arising chains of equations for 
(E, B, D) operator a valued functions ( the matrices of the finite dimensions 
for instance) coincides with the investigated before matrix Toda chain. Ex- 
plicit solutions for this chains of equations with the fixed ends reader can find 
in 0. Uniting these results it is possible to construct multi soliton solutions 
of the matrix 3-wave problem in three dimensions similar to way proposed 
in for construction of multi soliton solutions for matrix Devay-Stewartson 
equation. 

7 Outlook 

The concrete results of the present paper are concentrated in explicit formulae 
for discrete transformations for 3-wave problem of the section two and their 
generalization on the matrix case (section 6). 

But a no less important is the understanding how the method of the 
discrete transformation may be generalized on the case of multicomponent 
systems, connected with the semi simple algebras of the higher ranks r. 
From results of the present paper it is clear that in the case of arbitrary semi 
simple algebra there are r independent basis mutually commutative discrete 
transformations. In what connection are this commutative objects with the 
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main ingredients of the representation theory of the group is very interesting 
and intrigued question for further investigation. 

And the last comment. The chain with two fixed ends can not be con- 
sidered as the basis for some finite dimensional representation of the group 
of the discrete transformation, if it is at all possible to apply term group for 
it in this case. On the function at the end point of chain it is impossible 
to act by direct transformation at the right side and inverse on the left end. 
What is discrete transformation from the group theoretical point of view in 
this case? We at this time have no answer on this question. 
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9 Appendix I 

In this appendix we would like to show, how it is possible to construct two- 
dimensional integrable systems connected with A 2 algebra. 
Let us consider the following 3x3 polynomial matrix 



where are the polynomials of the degree k (with respect to parameter 
A and sign P means that coefficient function on the highest degree of the 
corresponding polynomial equal to unity. 

Let us define coefficient function of all polynomials are defined from the 
condition that between the columns of the matrix 



takes place the linear dependence in (rii + n 2 + + 3) points of the A plane, 
hi, h 2 Cartan elements of A 2 algebra and T{ = 4>i{t, A) + fi(x, A). fa are 
arbitrary rational functions with respect to argument A. 




(10) 



P = P(\) exp(n/ii + r 2 h 2 ) 
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The last condition is equivalent to the following system of linear equations 
for defining the coefficient function (we present it here for elements of the 
first row) 

Pnl + i(K) + c s exp(r| - 2r 1 s )F n 1 2 2 (A s ) + d s exp -(r| + rf)P£(X s ) = (11) 

8 = 1,2, («i + n 2 + n 3 + 3), r? = r(X s ) 

(ni + n 2 + n 3 + 3) exactly the number of coefficient function of polynomi- 
als of the first row. Thus (O) is the linear system of equation for their 
determination. 

Let us now determinant Det(P(X) = Det(P(X)). From ( [TCP it follows that 
it is the polynomial of (ni + n 2 + n 3 + 3) degree with unity coefficient before 
the highest term and from condition (|TT|) that it has zeroes in (ni+n 2 +ri3 + 3) 
points A s of the A plane. Thus 

Det(P(X)) = Det(P(X)) = n£L\ +n2+ri3+3) (A - X k ) (12) 

Now let us calculate the matrix PP~ l , where / means the differentiation 
with respect to one of two independent arguments of the problem x, t. From 
the definition of the inverse matrix it follows that matrix elements of this 
matrix are the following ones 

/fip-iN _ Det{P p -> P a + P a (i i+1 - fj) ^_^_ n n „s 
\rr )a,p ~ (ni+n2+n3+3) — , t - r 3 - U {16) 

This symbolical form means that the determinant matrix of numerator arises 
after exchanging of the (3 row of the P matrix on the a row of the matrix 
Pexp -(ri/ii + T 2 h 2 ). 

It is not difficult to understand that matrix PP~ l possess all the same 
singularities as functions r by themselves. 

Now let us illustrate situation on the example of three wave interaction, 
choosing T\ = X(cit + c 2 x), r 2 = X(dit + d 2 x) . Let us calculate in this case for 
instance (P t P~ 1 )n. In connection with ( pT3|) we numerator determinant have 

Pni+l + Pn x +lClX P n2 + P„ 2 (c 2 — Ci)A P nz — P n . A C 2 X 
D > 1 \ P ni Pnz+1 fn 3 

P-ni Pri2 Pns+l 
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It is obvious that between the columns of the matrix P the linear dependence 
takes place with the same coefficients and so numerator determinant has 
zeroes in the same points as determinant in enumerator. Computation the 
degrees of numerator shows that it is polynomial of the (ni+712+713+4) order 
and so considered matrix element is the linear function of the A parameter. 
From ( |i~2|) it follows that it equal exactly C\X. The same not conversion 
calculations show that matrix PtP~ l coincides with the u matrix from the 
introduction after identification 

p=(p 12 )2, Q = (p 13 )2, b = (p 21 )z, 

A D=(P 3l ):i, E=(P 32 )Z 

where values above are coefficients at the highest degree terms of the corre- 
sponding polynomial. These terms are known from the solution of the linear 
system (|ll|) and so we have explicit solution of the system (p. 
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In this appendix we would like to consider the simple example of soliton 
solution of 3-wave problem. We specially consider this simplest example 
in details to give the reader possibility to feel self-consistent of the whole 
construction of the present paper. 

Let in notations of the 5-th section £3 = 2, l\ — 0. Condition that vector 
I)2,o = is equivalent to the following system of equations 

A 3 = Af = A! = (14) 

The first of this equations leads uniquely to explicit form of initial D function 

D = 1 (l)/ 1 (2) + 2 (l)/ 2 (2), <j> x = <f/, f 2 = f (15) 

Using the initial conditions (f|) equation Af = may be rewritten conse- 
quently 
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(/1/2 - hh)Det I <P[ 0' 2 </i + &7a 

Keeping in mind that 0, / are the functions of the different arguments we 
conclude the last equations is equivalent to equality to zero of the two de- 
terminants of third order. The last condition in its turn can be rewritten as 
the system of equations 

0i = P01 + #2, 02 = «01 + ^02 

01=P0'l + #2, 02 = S0' 1 + t0 2 (16) 

0i" = p4>i + q&'L 02 = s 0'i + ^02 

From (PT6[ ) it follows immediately that (02 7^ c0i) p' = q' = s' = t' = and 
functions 0i )2 are the solutions of the first row of (pT6|)-th.e linear system of 
equation with the constant coefficients. Solution of this system is obvious 

</>i = ci exp An + c 2 exp A21 , 02 = c 3 exp An + c 4 exp A21 . 

From the equation EAf by the same way we obtain 

fx = di exp fi 12 + d 2 exp ^22 , h = <h exp Hvz + d A exp ji 22 , 

where c, d, A, /i arbitrary numerical parameters. 
The initial conditions 

—D 3 = D x + D 2 = BE = (b x exp/ii2 + b 2 exp fj, 2 2)(e 1 exp A x l + e 2 exp A 2 1) 



allow using ( |i~5| ) allow determine parameters b, e and find one relation con- 
nected parameters c, d, A, \i. Now let us calculate vector Q 2 ' using explicit 
expressions for D,B, E functions. The last two ones we present in the fol- 
lowing form E = p0i + q<p 2 , B = rf\ + 5/2- 



Q 



2.0 D 0i/i + 2 /2 



0'i 02/ V/2 /a 



D (A p\ D (4n 02 

p 2,o = V/2 qJ ^2,0 V r s 



D 1 fl fl\ jj I " 1 



h h) ~ V0i 02 
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Conditions of reality leads to other restriction on parameters involved. It is 
clear that two possibilities in the choice of parameters A and fi are possible 
A 2 = — A*, Ai = — Ai* and A 2 = — A 2 *. And the same limitations on parame- 
ters /i 2 - We do not present here explicit form for the other restrictions. This 
is pure algebraic manipulations. 
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